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A EEEL3MINAEY THEORETICAL STUDY OF AERODYNAMIC 
INSTAB ILITY OF A TW0-3LADE HELICOPTER ROTOR 
By Robert P. Coleman end Carl V. Stamp In 


SUMMARY 


A tlieory has lean developed in prelim: nary form which, seems 
capable of predicting the aerodynamic instability phenomena of a 
two -blade see-saw-type helicopter rotor. In particular, the theory 
indicates .the possibility of unstable Tibrations even with the 
chordwise center of mass at or ahead of the 25 -percent- chord 
position. 

The stability condition for oscillatory motion is expressed 
in terms of a small number of composite parameters that are evalu- 
ated from the moments of inertia, angle settings, and aerodynamic 
parameters of • a blade.. 

Computed stability results for different coning angle settings, 
center-of-mass positions; and control- system stiffnesses for 
one value of blade density and aspect ratio are presented in a 
chert. 

It 1 b found that, in addition to parameters analogciis to 
those occurring in wing-flutter theory, the present the ,\ry contains 
a parameter that represents' an. unstabilizing effect due tc the dif- 
ference between the moments, of inertia in flapping and- in rotation. 


INTRODUCTION 


TIlo present preliminary paper presents eon© numerical results 
of flutter calculations for. a two-blade helicopter. The results 
are published in preliminary form in order to male them more quickly 
available for possible application to the study of the phenomenon. 
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of "blade "weaving, " •which is apparently an aerodynamic instability 
and has bean called weaving from the appearance of the wavy path 
traced by the blade tips. The general method, which is an exten- 
sion of the method of reference (l), is to be treated in more 
detail in a later report. Derivations and discussions of the* method 
have consequently been largely emitted in the presont paper. Several 
effects are omitted which undoubtedly influence the quantitative 
results but are believed to bo unessential in a preliminary study 
of trends . 


METHOD OF ANALYSIS 


The computations are intended to apply to a helicopter having 
two rigidly connected blades sot at a coning angle and pitch sotting 
with respect to each other. The combined blades arc treatod as a 
singlo rigid body having three degrees of freedom in rotation about 
a fixed point at the hub. 

Tho air forcos aro obtained from a blade -element analysis 
using the wing flutter thoory of r of or once 2 but using tho theo- 
retical stoady state value for tho slope of tho lift curvo instead 
of tho complex function F + iG. The effect of a flexible control 
system is represented by a spring in the blade feathering degree 
of freedom. 


In tho air-force torms, the noncir culatory and the circulatory 
terms are considered separately. The noncirculatory tana are 
treated as apparent added mass and moment of inertia and combined 
with tho actual blade mass to obtain resultant blade, cantor- of- 
mass and moment -of- inertia parametors. 

The equations of motion are obtained by first writing tho 
Euler equations of motion for each blade, considered as a rigid 
body rotating about a fixed point. The air-farcc terms are expressed 
in terms of the angular-volocity compononts that occur in Euler's 
equations. 

The equations for tiro blades aro then combined in a way to 
represent a single rigid body.- The two blades aro Jointed at 
fixed angles with respect to each other and in such a way that the 
respective lines of center of mass intersect in a point at tho hub. 
Tho coordinate axes of Euler's equations are then transfo.ocod to 
axes rotating uniformly with tho moan blade motion. A spring- 
stiffness term is then inserted to represent tho effect of tao con- 
trol system.. 
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Vilen the various physical parameters have been combined into 
a smaller number of composite parameters to he defined in the fol- 
lowing section of the report, the equations of motion in matrix 
form become: • . , 


q + e d Ia, 


' ‘ 5 + H 3 


Q + % 


a 3, o 1 

Q -i 0 


1 k'° 

Cl 


0 0 rn 1 

L. QJ 


%» 0 0 

0 0 0 

0 0 0 


= o 


where 

e, P, r 

«s 




feathering, flapping, and lagging angpLes, respectively 


mean angular velocity of rotor 
and the other eytabcrls are, defined' in the section, Physical Parameters. 
If solutions are assumed of the form. 


a . IQfc 

. 0 ss A^e 


P - A p e ' 


XQb 


r-V® 

whqro. Ag, Ap, Ay are constants, the determlnartal equation for 
obtaining the values of X is 


/ 
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X 2 + E^X + I A + Kg' -x(l - I A J - % 0 « 

X(l - I B ) + H b X 2 + %X + I b 0 =. 0 

0 0 X 2 + H^X 

The equation shows that the y degree of freedom is independent 
of the others and need not he considered in the stability computa- 
tions. 

If the determinantal equation is expanded, a quartic in X is 
obtained of the form 

X^ + aX 3 + bX 2 + cX + d » 0 (l) 

where 

a = % + % 

b = I B + HgHjj + X A + Kq* + (l - I A ) (l - I B ) 
c = IjEq + ^I A + V>% + - I B )% + (l - x a )h b 

d = (i a + k^»)i b + HbH d 

for which the Important stability condition far oscillatory motion 
is 

c 2 - abc + a 2 d < 0 


-WV S + (% + %) (i-H^-IaI^Eb^* - (% + Hp) 2 (l-I B )V < 0 
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Kile can "be -written 


b 


' K e t ' s s^D 





± - *A 

go 

3 b 


1 + % 



( 2 ) 


The critical condition for divergence Is recognized, from equa- 
tion (i) a a 4s0j belt no further discussion of this type of 
instability will be given in the present paper. 


HEiSICAl FABAHEX3BS 
General Case 


The final stability conditions are engrossed in terms of 
certain composite parameters that are obtained frcaa ether basic 
parameters as fellows: Let A, B ; C denote the principal 

moments of inertia of a single blade (in a vacuum) about radial, 
choruvise, and perpendicular cues, respectively. The apparent 
added mass effect changes B to £<(1 + ft), where K is the 
ratio of the mass of a cylinder of air of a diameter equal to the 
chord of the LJnde to th-s mass of the blade, both eaten for equal 
length- along the span. A suitable average value applies to tapered 
blacs.es » . , 

The effect of apparent added maBS on A and C has been 
neglected, in the present computations. 

The following parameters involve the circulatory air force 
expressions for a single blade: 


« 
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■where r is the span-wise coordinate along a blade,, B is the tip 
radius, and p, b, a have the same meaning as in “wing flutter 
theory (reference 2), namely: 

p mass of air per unit of volume 

b ha .If chord of blade 

a chordwise distance of elastic axis (assumed to coincide -with 
center of mass) aft of midchord, divided by b. Leading 
edge is -Ij trailing edge is 1. 

Por the two blades rigidly connected so as to have a coning angle p Q 
and collective pitch setting 0 Q the corresponding parameters, with 
apparent mass included, are: 


A* = A cos 2 P q + B(l + k) sin 2 0 o sin 2 0 o + C cos 2 0 o sin 2 0 o 
B* = B(l + k) coq^o + C si n% rt \ 

v O i 


C* = A sin 2 p 0 + B(l + k) sin% 0 cos 2 P 0 + C cos 2 0 q cos 2 0 o J 
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H 1* = H 1 003 e o 003 P 0 ’- H 2 3±tt & o Bin ( cos 29 o 003 Po + 2 003 6 o Bin 

i 

Bin. fi Q sin p Q cos P Q + sin S £ 0 

Eg* = Eg cos e 0 ^cos 2© 0 cos - 2 sin 2 8 C sin ‘ 

p & ! 

= -2H-J. sin © Q cob © Q sin & Q cogfl 0 + cob Q q cob & 0 cos 2fl 0 i 

E^* * -2Eg sin 0 Q Bin 0 O cos 3 0 + cos^Pq 


then 


l B 


C 1 - B» %1 
A* + A* 


- c * - At _ J %* 
B r ” 3* 


V 


%“3? 


> 


^ “ 


V 

A* 


*e 


t _ 


A«n 2 


(6) 


■where Kg la the spring constant in the feathering degree of freedom. 


Special Cases 

Bor a homogeneous rectangular "blade the espressiona in equa- 
tions (3) "became 
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H-L ••= -2a KB 


H~ = 


in. = 


IT- = 2i 


-P$B 
2 b 


3=(| - a) k|a 

(M« 


(7) 


ileo 


B __ r£ 

■a' ^2 


(8) 


If the blade mass Is distributed la the plane of the span and chord 


A + 3 - C 


If P q 3 and other third-order angle tercas and certain. other com- 
binations of small terns are neglected, equations (4) and (5', 
become 


A* = A + Bp Q 2 
B r = £(1 + tt) - k£9 0 2 
C* = C . 3P 2 Bk0 


V = 


c. 


o 

<> 0 * * s e - ' 


7 - Wo - Wo + \ 


V " Ha(i - 


f' 59 + P 0 2 \ 


F( . ... gg , p . xr , A 5Pn- + 6 G "\ 

hi. -hhi- 0 P 0 + htj.{l - h t 

\ - / 


Hj= „ JB*£ q + F 5 (j. - 2;, 0 2 ) 


(9) 
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STABILITY CHART 


Stability conditions in terms of physical parameters can now 
he computed from equations (2) and (IQ) . Figure 1 is a plot of 

a/fc against {J for fixed values of the other parameters. 

The quantity j/k^/A may he recognized as the natural frequency 

of the control system and hlade in feathering motion at zero 
coning angle. 

The following constant values were assumed in the calculations: 

- « 1000 
A 

f oo 

K = 0.018 
% - 0.806 


Those values could he realized approximately, in a hlade having 
the following characteristics: 


Kadius, feet 23 ’- 75 

Chord, feet 1.5 

Weight, pounds 180 


The value 0=0 has, for simplicity, heen chosen in thl3 first 
calculation. An Inspection of the equations indicates that positive 
pitch angles 0 make the rotor more unstable. The effect of the 
chordwise center of mass upon the stability Is Indicated by the 
different curves. It is noted in particular that instability can 
occur even if the center of mass Is ahead of the 25-par cent- chord 
position. 


DISCUSSION OF EBSULTS 


The present theory appears to check, In a general way, the 
observations that instability of a two-blade see-saw rotor with a 
coning angle can occur even with the center of maso at or ahead 
of the 25-percent- chord position. 
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The theory shows that instability depends strongly upon coning 
angle and tc a lesser extent upon pitch sotting. A study of the. 
parameters employed in the analysis shows that and Eg choogs 

over a wide range with chants of coning angle and pitch setting, 
whereas Ig and Eg are relatively insensitive to angle changes. 

It appears after soma experience with this type of analysis that the 
parameters 1^ and Eg . are useful one 3 for designers to study in 
considering the effect of various design changes upon stability. 
Changes that increase the positive value of Ij^ improve the stability. 
Examples are forward movement of the chordwise center of mass, 
decrease of coning angle, and decrease of aspect ratio. The 
parameter Sn arises from a term representing aerodynamic damping 
m feathering oscillation. This type of tern is frequently ignored 
in studies of rotor-blade dynamics, but it seems to be important 
fer evaluating the effectiveness of increased control system 
stiffness in overcoming instability. If Eg is arbitrarily allowed 
to approach zero in this stability condition (2), this condition 
reduces to 

‘ I a)< 0 

which doeB not contain Kg : . It can be shown that the curves of 
figure 1 then became vertical straight linos and erroneously 
predict that stability is independent of rotor speed and control- 
system stiffness. 


Condition (2) shows that if a rotor is stable for 3g* = 0 
it will be stable for all values of Eg*. These points, corre- 
sponding to Kg' = 0, therefore represent important suff icier. t 

conditions fer stability. If Kg* is allowed to approach zero 
in condition (2), the sufficient condition for stability becomes 


l B 


1 - IA 

h d/% 


- 1 + % - V 


<0 


The part of this expression that is sensitive to angle changes is 
the parameter 
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which can. be expressed in the form 

A 1 + B* - C* + Hi* 

P — 

' V 


For stability the parameter P should not Toe larger than a certain 
positive quantity determined mainly by Hp. Changes in a ship that 
decrease the positive magnitude of P will improve the stability. 
If P is evaluated in terms of angles from equations ( 10), it 
becomes, with certain small terms neglected 

(1 + 2a) + 2 ^ + 2 e p E 

P - * ± g k 

[fl \ E 

-a 3^ - -46<A>7 

_ — 


HOW TO PREVENT 2SHEABIIJTI 


Insofar as the present theory has included the correct degrees 
of freedom and types of forces to account for the phenomenon of 
instability it can be used to indicate the effect of proposed 
remedies, such as the increase of control- system stiffness, forward' 
displacement of chordwlse center of mass, decrease of coning angle, 
or decreaee of aspect ratio. 

The general result of the analysis ia that a see -saw rotor with 
a coning angle is more unstable than an airplane wing having 
corresponding parameters. The additional unstabilizing effect is 
associated with the difference in moments of inertia in flapping 
and in rotation, A rotor is most stable when the moment of inertia 
in rotation C* is larger than the moment of inertia in flapping B f . 
The main effect of increased coning angle appears to be to decrease C f 
Without making a compensating change in B* . Conf igurations that 
tend to confine the mass distribution to the plane of rotation ere 
therefore desirable for stability, Configurations that tend to 
spread out the mass in the plane of flapping are undesirable ter 
stability. 
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The effect of mechancial damping lias not teen explicitly 
treated, but it is expected to have a stabilising effect by general 
analogy with wing flutter theory. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Ya. 
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